Topological Classification of Crystalline Insulators with Point Group Symmetry 
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We show that in crystalline insulators point group symmetry alone gives rise to a topological 
classification based on the quantization of electric polarization. Using C3 rotational symmetry as 
an example, we first prove that the polarization is quantized and can only take three inequivalent 
values. Therefore, a Z 3 topological classification exists. A concrete tight-binding model is derived 
to demonstrate the Z3 topological phase transition. Using first-principles calculations, we identify 
graphene on BN substrate as a possible candidate to realize the Z3 topological states. To complete 
our analysis we extend the classification of band structures to all 17 two-dimensional space groups. 
This work will contribute to a complete theory of symmetry conserved topological phases and also 
elucidate topological properties of graphene like systems. 
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Since the celebrated discovery of the quantum Hall 
effect [1 , topological classification of electronic states 
has emerged as a powerful concept in condensed mat- 
ter physics. The quantum Hall insulators are distin- 
guished from ordinary insulators by a topological index, 
the TKNN number, which gives the quantized Hall con- 
ductance [2j [3]. For a long time, the TKNN number 
was thought to be the only topological index describ- 
ing non-degenerate electronic ground states. Recently, it 
was realized that in crystalline insulators new topological 
indices can be defined in the presence of discrete symme- 
tries, which has led to the identification of a slew of new 
topological states. For example, the quantum spin Hall 
insulators are characterized by a nontrivial Z2 index [4], 
which is protected by time-reversal symmetry [5j[6]. Sim- 
ilarly, magnetic translation symmetry can also give rise 
to a Z2 classification in antiferromagnetic insulators [7]. 
Another interesting proposal is the so-called topological 
crystalline insulators, in which a Z2 index can be defined 
and is protected by both time-reversal and certain point 
group symmetries [3J[9]. 

In this Letter we show that in crystalline insulators 
point group symmetry alone can give rise to a new topo- 
logical classification based on the quantization of electric 
polarization. Our idea is inspired by a beautiful result 
due to Zak [TO] , i.e., in one-dimensional (ID) systems 
with inversion symmetry, the Berry phase of the Bloch 
bands can be either or tt [11]. This quantization of 
the Berry phase naturally leads to a Z2 classification in 
ID. In higher dimensions, we find that the role of Zak's 
phase is replaced by a closely related quantity, the electric 
polarization [12] [13] , which is quantized in the presence 
of point group symmetry. The generalization to higher 
dimensions is expected to display a richer spectrum of 
possibilities because of the enlarged symmetry class com- 
pared to ID. 

For the sake of definiteness, 2D crystals with C3 rota- 



tional symmetry are used as an example in the following 
discussion. We first prove, from general symmetry ar- 
gument, that the polarization is quantized and can only 
take three inequivalent values. Therefore, a Z3 topologi- 
cal classification exists. We then proceed by constructing 
a concrete tight-binding model to demonstrate the exis- 
tence of various Z 3 topological states and the topological 
phase transition in the presence of C3 symmetry. Using 
first-principles calculations, we identify graphene on BN 
substrate [T414T7] as a possible candidate to realize the 
Z3 topological states. The misconception that C3 sym- 
metry enforces zero polarization is also clarified. Finally, 
we extend our analysis to all space groups in 2D. Since 
a many-body formulation of the electric polarization al- 
ready exists p~8j [19] , the topological classification scheme 
discussed here can be readily applied to interacting sys- 
tems. 

General symmetry analysis. — Theoretically, it has 
been well established that the polarization in a crys- 
tal is formally defined only modulo a quantum uncer- 
tainty [12] [13], i.e., the polarization vectors P given in 
the following expression are all equivalent to each other: 

P-P + e^n^, (1) 

i 

where e is the electron charge, rii are integers, and di are 
the primitive lattice vectors. This is a direct consequence 
of the translation symmetry of the lattice. At the same 
time, P has to also satisfy the point group symmetry of 
the crystal, which places great restriction on the allowed 
values of P, thereby quantizing it. For example, in ID the 
only point group symmetry is inversion, which requires 
that P and — P must be equivalent. This leads to two 
possibilities of P (0 and ea/2) and Zak's phase (0 and 
7r) [10], and a Z 2 classification. 

Let us now consider a 2D lattice with C3 rotational 
symmetry. We denote the two primitive lattice vectors 
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by a and b. The polarization vector can be expressed as 
P = e(aa + /3b). Consider the in-plane rotation R by 
an angle of 27r/3. R operates on the lattice vectors such 
that 



R(a) = b — a , R(b) = —a , 
and P transforms according to 

R(P) = P' = e(ab- (a + 0)a) 



(2) 



(3) 



However, the transformed polarization vector P' must be 
equivalent to the starting value P according to Eq. 0, 
i.e., 



P' = P + e(ma + nb) 



(4) 



Combining the above equations of P' yields possible val- 
ues for a and f3: 



a = 



2n 



(5) 



where n, m and j are integers. The polarization values 
allowed by Cs symmetry thus fall into 3 inequivalent sets 
described by 0, (a + 6)/3, 2 (a + 6)/3, and corresponding 
series of equivalent values given by Eq. 0. 

This result has two important consequences. First, it 
shows that while zero is an allowed value, Cs symmetry 
does not preclude nonzero polarization values as is some- 
times assumed in the literature. This misconception is a 
result of the classical intuition that like forces, one can 
cancel polarization vectors in opposing directions. While 
this may apply for dipole vectors in finite, discrete charge 
distributions, it is not true for bulk crystalline insula- 
tors. Secondly, and more importantly, the above result 
suggests that 2D crystalline insulators with Cs symme- 
try can be distinguished by the three possible values of 
P, giving rise to a Zs classification. Following the same 
procedure, one can verify that Cq rotational symmetry, 
e.g., as seen in graphene, allows only zero polarization in 
insulators. 

Zs Topological phase transition. — We now provide a 
concrete example of Zs topological states by construct- 
ing a tight-binding model. Let us consider a honeycomb 
lattice with two orbitals, denoted by <j) a and <^g, on each 
site. We assume that there is a nearest-neighbor hopping 
t between orbitals from the same type, and an on-site 
inter-orbital coupling t' between the a and f3 orbitals. 
To lower the symmetry from Cq to C3, the A and B sub- 
lattices are assumed to have opposite site energies, ±A a 
and d=A/3, for each type of orbitals. The tight-binding 
Hamiltonian is given by 



(ij) 



(6) 




FIG. 1. (color online) Topological phase diagram of the 
tight-binding model described by the Hamiltonian |7|. The 
solid (blue) curve is the phase boundary given by A a = t' 2 . 
In each region, the polarization is given by Pi = — (a + &)/3, 
P11 = 0, and Pin = (a + b)/3. The dashed line indicates 
a representative path in the parameter space that undergoes 
the topological phase transition twice. 



where & = ±1 for A and B sublattices. The correspond- 
ing Bloch Hamiltonian is obtained by Fourier transform: 



H(k) = 
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where the structural factor is given by 14 = 1 + 
2 cos(v / 2/2^) cos(k y /2) - 2i s'm(V3k x /2) cos(k y /2). 

Next we show that at half-filling the system can be 
tuned to three topologically distinctive states by chang- 
ing the parameters without breaking the Cs symmetry. 
We first identify the gap closing point in the parame- 
ter space because P cannot change unless the band gap 
closes. The band energies at the K point is given by 



E = ±(A a + Ap)±J(A a -Ap)* + 4t> 2 . (8) 



Obviously, when t' = A a Ap the band gap will close, 
which signals a topological phase transition characterized 
by a sudden change of P. To confirm this statement, we 
calculate the polarization according to the Berry phase 
formula [T2j|13] at three representative points in the A a — 
A p plane, and indeed find three inequivalent values of P 
predicted by Eq. ([5|. The resulting topological phase 
diagram is shown in Fig. [l] 

Material realization. — We now present a realistic ma- 
terial system, graphene on /i-BN substrate [T4HT7] , in 
which various Zs topological states can be realized. 
Specifically, we consider four cases: fo-BN, AA stacking 
of graphene on fo-BN, and AB stacking of graphene on 
/i-BN with either boron or nitrogen atoms sitting directly 
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FIG. 2. (Color online) (a),(b) Schematic illustration of po- 
larization in systems with C3 symmetry, (a) /i-BN lattice 
with 3 equivalent unit cells, (b) Possible values of polariza- 
tion. Each set of topologically-equivalent values is marked by 
a specific symbol, (c)-(f) displays calculated values of polar- 
ization for different structures and the corresponding topo- 
logical classes, (c) Structure 1: /i-BN monolayer; (d) Struc- 
ture 2: A A stacking of graphene monolayer on /i-BN; (e) AB 
stacking of graphene monolayer on texth-BN with a carbon 
sublattice lying above the boron sublattice and (f) Structure 
4: AB stacking with a carbon sublattice lying above the nitro- 
gen sublattice. The unit cells and lattice vectors a and b are 
marked in black, whereas the polarization vector is marked in 
red. Colors used to represent various atoms are listed. 

underneath the carbon atom. Although graphene has Cq 
rotational symmetry, placing graphene on h-BN reduces 
the symmetry from Cq to C3, allowing non-zero polariza- 
tion values to emerge. 

We have performed first-principles calculations to ob- 
tain the polarization in these structures and compare 
with our symmetry analysis. While the carbon-carbon 
(C-C) bond length in graphene is known to be 1.42 A, 
/i-BN has a B-N bond length of 1.46 A [20]; hence the 
mismatch is very small ~ 2%. Our initial structures com- 
prised, as shown in Figj2j of bulk monolayer /i-BN and 
various stacking of graphene monolayer on it. We lattice 
matched /i-BN to graphene by taking an initial lattice 
constant of 1.417A and constructed super cells represent- 
ing these materials in bulk (see Fig{2|. The out-of-plane 
distance between graphene and /i-BN was taken to be 
3.2 A [14 . The vacuum along the z-direction was taken 
to be 30A approximately. We used the Vienna Ab-initio 



TABLE I. In-plane polarization values for the structures 
shown in Fig[5] specified in terms of lattice vectors a and 
b. We tabulate the values of the ionic part Pi on and the elec- 
tronic part Peie, followed by the total value of polarization 
P. All values are in terms of electronic charge x lattice vector 
such that P = e(aa + /3b). 
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Simulation Package (VASP) [21] which uses a plane- wave 
basis set and employed PAW pseudopotentials [22] , along 
with a local approximation to the exchange-correlation 
potential (LDA) [23 . We set an energy cut off of 400 
eV and a /c-mesh of 4x4x1 for the self-consistent run and 
36x36x1 for the subsequent polarization calculation. 

The total polarization can be divided into the elec- 
tronic and ionic parts. While the electronic part of the 
polarization was calculated using the Berry phase for- 
mula [12] [13] , the ionic part was obtained from the dipole 
density distribution. Polarization values thus calculated 
for bulk structures shown in Figj2j are listed in Table [I] 
in corresponding order. We list a and j3 for every P, such 
that P = e(aa + /3b), where a and b are the correspond- 
ing lattice vectors of the system. In complete agreement 
with our symmetry analysis, all obtained values fell into 
and exhausted the 3 topological categories predicted by 
Eq. (§. 

Here we comment on the zero polarization found in h- 
BN [24 , 25 . The electronic structure of /i-BN can be well 
described by taking the ionic limit, i.e., the boron atom 
loses all its valence electrons to nitrogen. The nominal 
charge at the boron and nitrogen sites are therefore +3 
and —3, respectively. (One can verify that in the ionic 
limit, the ionic and electronic contributions agree with 
the DFT calculation.) This guarantees that the total 
polarization of /i-BN will always be an integer multiple 
of a + 6, which is equivalent to zero. If we make an 
artificial crystal by replacing boron and nitrogen with 
beryllium and oxygen, then the total polarization will not 
vanish. So the zero polarization of /i-BN is not because 
C3 rotational symmetry forbids other values, instead it 
is due to simple electron counting. 

Polarization in 2D space groups. — Similar analysis can 
be carried out for all 17 two-dimensional space groups. 
The allowed values of polarization derived from symme- 
try constraints are shown in Fig. [3] 

Several remarks are in order. For symmetry group pi, 
the allowed value of P forms a continuum and fills up 
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FIG. 3. (Color online) Allowed values of polarization for all 
17 two-dimensional space groups. The corresponding Wigner 
Seitz cells for various lattices are displayed in blue. Allowed 
points and lines of points are marked in magenta. Allowed 
values of polarization are vectors joining the origin of the 
Wigner Seitz cells to the points marked. 



the entire Wigner Seitz cell. This is consistent with the 
fact the only point group symmetry operation of p\ is the 
identity operator and there is no topological classification 
protected by p\ . We also note that for several symmetry 
groups the allowed value of P forms a line in one direc- 
tion but is still quantized in the other direction. In this 
situation, a topological classification can be obtained by 
taking the discrete component of P. 

Summary — We have presented polarization as a topo- 
logical index protected by point group symmetry for crys- 
talline insulators. In particular, we derive 3 distinct 
topological classes for Cs rotational symmetric systems. 
Topological phase transitions between these classes are 
shown to be accompanied by closing of the band gap. 
We complete the analysis by deriving polarization values 
for all 17 2D space groups which gives rich possibilities 
of topologically distinct classes. The topological clas- 
sification scheme discussed here can be readily applied 
to interacting systems via the many-body formulation of 
polarization [THJ [19] . 
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Note added. — Upon the completion of our work, we 
noticed the preprint [26 on similar topics. 
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